Multi-Phonon 7- Vibrational Bands and the Triaxial Projected Shell Model 
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We present a fully quantum-mechanical, microscopic, unified treatment of ground-state band and 
multi-phonon 7-vibrational bands using shell model diagonalization with the triaxial projected shell 
model. The results agree very well with data on the g- and 7-band spectra in 
as with recently measured A'^ 2-phonon 7-bandhead energies in I'^^Er and ^"^^Er 
7-excitation energies are predicted. 
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The atomic nucleus is a many-body system with pro- 
nounced shell effects that can have intrinsic deformation. 
In addition, it can, according to the semi-classical col- 
lective model, undergo dynamical oscillations around the 
equilibrium shape, resulting in various low- lying collec- 
tive excitations. Ellipsoidal oscillation of the shape is 
commonly termed a 7- vibration [Q. 

Thanks to advances in high-resolution 7-ray detec- 
tors, high quality measurements not only of high-spin 
states but also of low-spin states are now commonly 
available. As a consequence long-sought multi-phonon 
7-vibrational states have been discovered in a series of ex- 
periments over the last decade |^|^]. However, the status 
of unified theoretical descriptions for ground-state band 
(g-band) and multi-phonon 7-vibrational bands (7-band) 
is not so satisfactory. In the present work, we attempt a 
consistent description of these low-lying bands using an 
approach based on the Projected Shell Model (PSM) 

In its original form, the PSM uses an axially deformed 
basis. The shell model diagonalization is carried out 
within the space spanned by the angular momentum pro- 
jected quasiparticle (qp) vacuum, 2- and 4-qp states. In 
this sense, the PSM is a Tamm-Dancoff approach and 
one expects that the collectivity of low-lying states may 
be strongly affected by mixing many 2- and 4-qp states. 
Indeed, a multi-qp admixture can cause significant effects 
in band crossing regions |^ . 

However, in the low-spin region before any band cross- 
ings (/ < 10), the admixture is very weak and the cal- 
culated g-band always exhibits the characteristics of an 
axially symmetric rotor. For example, the ordinary PSM 
fails to describe the steep increase of moment of iner- 
tia at low spins in transitional nuclei |^ . Quite recently, 
the restriction to an axially deformed basis in the PSM 
was removed by two of the present authors (JAS and 
KH). It was shown that the observed steep increase of 
moment of inertia for transitional nuclei can be well de- 
scribed if one introduces triaxiality in the deformed basis 
and performs 3-dimensional angular momentum projec- 
tion H. This approach is called the Triaxial Projected 



Shell Model (TPSM). 

Another important issue is whether the PSM can de- 
scribe bands built on collective vibrational states. The 
usual treatment of the 7-band based on the Tamm- 
Dancoff or on the Random Phase Approximation assumes 
different coupling constants for the /i = and — ±2 
parts of the QQ-force, with the former related to the 
mean-field deformation and the latter adjusted to the 7- 
bandhead energy. In the PSM, as in the ordinary shell 
model, such an adjustment is not permitted because the 
Hamiltonian must be rotation-invariant and thus these 
two coupling constants must be equal: One cannot sim- 
ply fit the theoretical 7-bandhead to the experimental 
one by modifying the QQ-force in that manner. 

On the other hand, one might hope that inclusion of 
many 2-qp states could introduce a collective contribu- 
tion that would produce the desired low-lying 7-state. 
But such attempts have failed. Because of the large pair- 
ing gap, the energy of the lowest 2-qp state is above 1.5 
MeV and is much higher than the actual 7-bandhead 
energy, which typically lies between 0.5 and 1 MeV in 
rare-earth nuclei. The QQ-force is too weak to lower 
the theoretical 7-band energy by such a large amount in 
a limited basis. Calculations including about one thou- 
sand 2- and 4-qp states do not lead to low-lying excited 
states that look like the experimental 7-band ||l^. One 
therefore has to conclude that it is not practical to de- 
scribe the 7-vibrational state in terms of multi-qp states 
in the framework of the axial PSM |pT[| . 

In the present paper, the TPSM extension of the PSM 
and the computer code developed in are used to study 
multi-phonon 7-bands. (Although the present theory is 
not based on a vibrational phonon excitation mechanism 
as in other models [|l2| , we shall use the conventional vi- 
brational terminology in our discussion.) We shall show 
the following: (1) For well deformed nuclei, introduc- 
tion of triaxiality in the basis does not destroy the good 
agreement for the g-bands obtained previously in the ax- 
ial PSM calculations (for example, those presented in 
Ref. B). (2) However, it produces new excited states 
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(7-bands) at the correct energies that do not occur in 
the axial PSM. (3) For transitional nuclei, use of a ba- 
sis of fixed triaxiality improves the g-band moments of 
inertia, as already shown in ||^, and at the same time 
produces realistic 7-bands. (4) By a single diagonaliza- 
tion of the Hamiltonian (with the same parameters in the 
deformed basis), we obtain not only the g- and 7-band, 
but also higher excited bands that can be identified as 
the multi-phonon 7-bands and these compare very well 
with recently measured 4"*" 2-phonon 7-bands. (5) Fi- 
nally, we make predictions for the 2- and 3-phonon 7- 
band (referred to as 27- and 37-band hereafter) of those 
Er isotopes treated here for which no measurement has 
yet been reported. 

Since an extensive review of the PSM exists (see Ref. 
1^ and references cited therein), we shall describe the 
model only briefly. The PSM (TPSM) closely follows 
the shell model philosophy and is, in fact, a shell model 
truncated in a deformed basis. One uses a Nilsson poten- 
tial having axial (triaxial) deformation to generate the 
deformed single-particle states. The Nilsson spin-orbit 
force parameters k, and ji are essential in reproducing 
correct shell fillings. For rare-earth nuclei, we use the 
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early compilation of Nilsson et al |13| without modifica- 
tion. For the axial deformation parameter e in the Nils- 
son model, we take the values given in Ref. Thus, 
for the TPSM, the triaxial deformation e' is the single 
adjustable parameter. The static pairing correlations are 
treated by the usual BCS approximation to establish the 
Nilsson-|-BCS basis. The 3-dimensional angular momen- 
tum projection is then carried out on the Nilsson-f-BCS 
qp-states to obtain the many-body basis, and the Hamil- 
tonian is diagonalized in this projected basis. 

In the present work, we consider only low-spin states 
where no band crossing with any multi-qp band occurs 
in the yrast region. Thus, the many-body basis may be 
restricted to the projected triaxial qp vacuum state: 
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where | > represents the triaxial qp vacuum state. 
This is the simplest possible configuration space for an 
even-even nucleus. Note that only one state is possible 
for spin J = (the ground state). Thus, multi-qp com- 
ponents have to be taken into account if one wants to 
describe 7 = excited states (see further discussion be- 
low). The diagonalization is performed over a chain of 
Er isotopes up to spin / = 10. 

As in the usual PSM calculations, we use the Hamil- 
tonian [|| 

H = H,- \xY.QlQ, - ^^^P^P -GqY. (2) 



so that the corresponding Nilsson Hamiltonian (with tri- 
axiality) is given by 



Here ifg is the spherical single-particle Hamiltonian, 
which contains a proper spin-orbit force as mentioned 
before, while the interaction strengths are taken as fol- 
lows. The QQ-force strength x is adjusted such that the 
physical quadrupole deformation e is obtained as a re- 
sult of the self-consistent mean-field (HFB) calculation 
The monopole pairing strength Gm is of the stan- 
dard form Gm = [21.24 =F 13.86(7V - Z)/A] /A, with "-" 
for neutrons and for protons, which approximately 
reproduces the observed odd-even mass differences in 
this mass region. This choice of Gm is appropriate for 
the single-particle space employed in the PSM, where 
three major shells are used for each type of nucleons 
{N — 4,5,6 for neutrons and N = 3,4,5 for protons). 
The quadrupole pairing strength Gq is assumed to be 
proportional to Gm, the proportionality constant being 
fixed as usual to be in the range 0.16 - 0.18. These inter- 
action strengths are consistent with those used previously 
for the same mass region [||J^. 

Let us first consider a well-deformed nucleus ^^*Er, 
which is generally considered to be axially symmetric. 
In fact, previous (axial) PSM calculation for this nucleus 
gave an excellent description of the yrast band up to a 
very high spin ||]. Fig. la shows the calculated energies 
as functions of the triaxiality parameter e' for angular 
momenta up to J = 10. In addition to the usual g-band 
with spins / = 0, 2, 4, • • •, a new set of rotational states 
with spins / = 2, 3, 4, • • • appears. This figure looks sim- 



ilar to the one shown by Davydov and Filippov |15|, but 
now obtained in terms of a fully microscopic theory. Un- 
like the irrotational flow model, the PSM spectrum de- 
pends not only on the deformation parameters but also 
on the shell filling of the nucleus in question. We see 
that, for the g-band of ^^^Er, the energies as functions 
of triaxiality are nearly flat and their values remain close 
to those at zero triaxiality. Thus, the triaxial basis has 
no significant effect on the g-band for a well-deformed 
nucleus and does not destroy the good g-band result ob- 
tained with an axially deformed basis. 

However, it has a drastic effect on new excited bands 
(second and higher excited bands are not shown in the 
figure) . Their excitation energies are indeed very high for 
axial symmetry, but come down quickly as the triaxiality 
in the basis increases. At e' — 0.13, the first excited band 
reproduces the observed 7-band in ^^*Er (while preserv- 
ing the good g-band agreement). It should be noted that 
the excited bands studied in this paper are obtained by 
introducing 7-degree of freedom in the basis (quasipar- 
ticle vacuum). They are collective excitations, but not 
quasiparticle excitations. We may thus identify the first 
excited band as the 7-band, the second excited band as 
the 27-band, the third excited band as the 37-band, etc. 

The above results can be understood by studying the 
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JC-mixing coefficients for each, projected if-state (see Eq. 
(|l|)) in the total wavefunctions. It is found that for this 
well-deformed, axially symmetric nucleus, if-mixing is 
negligibly small. States in the g-band are essentially the 
projected K = state for any e'. That is why the basis 
triaxiality does not destroy the result obtained with an 
axially deformed basis. The excited bands are also built 
by rather pure projected -fC-states. For example, the first 
excited band with the bandhead spin / = 2 is manily 
the projected K = 2 state and the second excited band 
with the bandhead spin / = 4 is the projected K = 4 
state. A small amount of if -mixing can be seen only for 
states with higher total spin if triaxiality in the basis is 
sufficientlly large. 

Fig. lb illustrates another example, the transitional 
nucleus ^^^Er. We see that the energies for the g-band 
are no longer constant, but clearly vary as functions of 
triaxiality. This feature is expected for a 7-soft nucleus. 
For the excited bands, triaxiality in the basis has a sim- 
ilar effect as we have seen for well-deformed nucleus: it 
drastically lowers their energies to those of the observed 
7-band. 

A rather different picture of ii'-mixing is observed for 
this 7-soft nucleus. The states are no longer pure pro- 
jected iiT-states, but highly mixed. For example, the two 
1 = 2 states (the one in the g-band and the other one be- 
ing the bandhead of the first excited band (the 7-band) 
are mixed from the projected K = and K = 2 states. 
At e' = 0.13, the 1 = 2 state of the g-band is contam- 
inated by the projected K = 2 state with a weight of 
about 1/4, and the 1 = 2 state of the first excited band 
contains the projected K = state with a weight of about 
1/4. Stronger JsT-mixing is seen for states with higher to- 
tal spin and larger basis triaxiality. 

Fig. 2 presents results for a chain of Er isotopes with 
neutron numbers from = 88 to 102. This covers both 
transitional {N w 90) and well-deformed {N « 98) nu- 
clei. The theoretical results are compared with available 
data for both g- and 7-band up to / = 10. The axial 
and triaxial deformation parameters used in the present 
calculations are listed in Table I. The triaxial parameter 
e' = 0.13 giving the correct position of the 7-band for 
^^^Er and ^^^Er corresponds to 7 = 25.5° in terms of 
the usual gamma parameter, if one uses as a very rough 
estimate 7 ~ tan~-'^(e'/e)- 

A microscopic description of transitional nuclei has al- 
ways been challenging. The nuclei discussed here with 
neutron number around 90 have g-bands that are quasi- 
rotational but with considerable vibrational character. 
The ground-state energy surface of a transitional nucleus 
was shown to have a shallow minimum at a finite 7- 
deformation in HFB calculations 1 16|. It has been demon- 
strated that such a shallow minimum becomes a promi- 
nent minimum when projected onto spin 7 = JT^ . The 
necessity of introducing triaxiality in the PSM basis to 
describe the observed g-band moment of inertia in tran- 



sitional nuclei was demonstrated in Ref . [gj . We now see 
that, with the same triaxiality, the first excited TPSM 
band reproduces also the observed 7-band. By adjusting 
a single parameter e' in the TPSM, the spectra of both 
the g- and 7-band are described simultaneously and con- 
sistently by the Hamiltonian (|^) diagonalized within the 
Hilbert space (^. 

Next, let us turn to a discussion of multi-phonon 7- 
bands. In Fig. 3, we plot all the states for spins / < 
10 obtained after diagonalization within our projected 
triaxial basis for two nuclei in which a 27-band has been 
reported. For ^^*Er, the second excited theoretical band 
agrees beautifully with the new 4+ 27-band reported in 
Ref. §]. For ^^^^Er, the observed 4+ 27-bandhead |] 
is also well reproduced. Since our theory agrees very 
well with the g-band and the (1-phonon) 7-band observed 
in these nuclei, the present results support strongly the 
interpretation of these data as 27-bands. 

To our knowledge, no 37-band has yet been seen ex- 
perimentally. According to our calculations, they should 
appear between 3 and 3.6 MeV. In Table I, we list the the- 
oretical values for the 2+ 7-, 4+ 27- and 6+ 37-bandhead 
energies. As Table I shows, the predicted 7- vibrational 
spectra are quite anharmonic. Anharmonic 7-vibrations 
have been discussed by several authors [p^|~pO|. This 
anharmonicity is a straightforward consequence of the 
present microscopic theory. This may be contrasted with 
earlier models that found it necessary to introduce ex- 
plicit anharmonicities to reproduce the 7-band spacings 

Finally, we mention briefly the 0+ excited states. Un- 
like the usual collective models based on phonon exci- 
tations, a O"*" collective excited state does not exist in 
the present calculation. Excited 0+ states can occur if 
we include multi-qp states on top of the present vacuum 
configuration. However, since the states constructed in 
this way are mainly qp in character, the collectivity of 
such a 0+ excited state is generally expected to be much 
weaker than that of a 2-phonon 7-state, which should 
have a large E2-decay probability to a 1-phonon 7-state. 
Furthermore, such states should depend strongly on the 
shell fillings. Therefore, the nature of the 4+ 2-phonon 
excited state is kinematical while the 0"*" 2-phonon ex- 
cited state is dynamical. There has been one experiment 
reporting a 0"*" excited state in ^^^Er ||^; the measured 
B(E2:0''"^2+) is enhanced, suggesting that this 0+ ex- 
cited state is of 2-phonon nature. At present, this is 
a single observed example of the 0+ 2-phonon excited 
state. Whether this observation can be reproduced by 
the TPSM with inclusion of qp states remains to be seen. 

To summarize, we have applied the Triaxial Projected 
Shell Model to some Er isotopes to investigate multi- 
phonon 7 vibrational bands. The shell model diagonal- 
ization is not carried out in a spherical basis as for a 
conventional shell model, but in a deformed basis with 
triaxiality. It is found that this simultaneously improves 
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the description of the g-bands in transitional nuclei and 
leads to a consistent description of multi-phonon 7-bands 
in both transitional and well-deformed nuclei. The newly 
observed 4+ 27-bands are reproduced by the same cal- 
culation, thus supporting their experimental assignment, 
and the bandhead energies of as yet unobserved 6'^ 87- 
bands are predicted. 

Thus, our unified view of the g- and multi-phonon 7- 
bands agrees surprisingly well with the existing data, 
even though we have used the simplest possible config- 
uration space. The origin of the 7-bands discussed in 
the present paper is kinematical rather than dynami- 
cal, indicating a microscopic connection between the 7- 
excited states and the nuclear ground state properties. 
We are presently investigating various intra- and inter- 
band B(E2)-values to test the theory further. These re- 
sults will be discussed in terms of X-mixing and reported 
in a longer paper. 

Dr. Kenji Hara^ worked on the present paper until his 
last day. This Letter is dedicated to the memory of his 
lifetime contributions to the Projected Shell Model. This 
work was supported in part by Conacyt (Mexico) . 
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TABLE I. Axial and triaxial quadrupole deformation pa- 
rameters e and e employed in the TPSM calculation for Er 
isotopes. 7527(4^) and £3^(6^) are predicted energies for the 
7-, 27- and 37-bandhead in units of MeV. Note their anhar- 
monicity. 
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